Detection of A^-particle entanglement with generalized Bell inequalities 
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We show that the generalized Bell-type inequality, explicitly involving rotational symmetry of 
physical laws, is very efficient in distinguishing between true A'^ particle quantum correlations and 
correlations involving less particles. This applies to various types of generalized partial separabilities. 
We also give a rigorous proof that the new Bell inequalities are maximally violated by the GHZ 
states, and find a very handy description of the A^-qubit correlation function. 
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I. INTRODUCTION 

The correlations of quantum systems containing more 
then two particles are an extensive field of theoretical 
and experimental study. Before the trailblazing paper of 
Greenberger, Horne and Zeilinger |1J in 1986 Svetlichny 
derived a Bell- type inequality [3, violated by quantum 
mechanical predictions for pure three-qubit states, with 
the property that it can distinguish between three-qubit 
correlations and two-qubit ones. A generalization of the 
results by SvetHchny was given in e.g. In Q one 

can find an inequality for N qubits, which can be used 
to classify the iV-body correlations in the following way. 
If any quantum state violates the inequality by a factor 
larger than 22(^^2)^ then the state contains true A^-body 
correlations. The factor that distinguishes maximal vio- 
lation of the inequality, which occurs for a GHZ state, and 
the maximal possible violation by a state with one qubit 
completely unentangled with other — 1 ones is V2. In- 
terestingly in Ref. 5] it was shown that a suitable choice 
of measurement settings allows to increase this factor, for 
three qubits, up to two (this is important, because the vi- 
olation factor is directly linked with the resistance, with 
respect to noise admixtures, of the non-classical correla- 
tions) . Here we investigate the Bell inequalities of a type 
introduced in Q , which include the assumption that 
physical laws must have a rotationally invariant form. 
They are dependent on the measured values of the corre- 
lation tensor for the A^-qubit processes under investiga- 
tion. What is to be measured are 2^ components of the 
tensor, thus the same number of measurement settings is 
involved as for the standard Bell inequalities. The fac- 
tor of violation of such an inequality for a GHZ state 
is two times higher than for any biseparable state. The 
results are very easily generalizable to more complicated 
situations. We also give a rigorous proof that the Bell in- 
equalities are maximally violated by the GHZ states, and 
find a very handy description of the A'^-qubit correlation 
function. 

As it was shown in Ref. 7J the generalized Bell in- 
equalities can be utilized to reveal to what extent realistic 
models can reproduce correlation functions E{ai, a^r). 



which have a rotationally invariant form. It is assumed, 
that correlation functions are linear in vectors a^, which 
specify the dichotomic local qubit observables (i.e., di - a, 
where a stands for the three Pauli operators) . Note that 
such is the case for the quantum correlation functions. 

Let us present the basic elements of our reasoning. On 
one hand, in the case of local and realistic theories, the 
correlation function for the qubits must be expressible in 
the following form: 



Envidiiai), . 



N 

,dN{aN))= / dAp(A)]^/i(ai), 



(1) 



where p(A) is a certain distribution function of some (hid- 
den) parameters A, and li is a function, that predeter- 
mines the values of results of experiments that can be 
performed on the given local system. Its allowed values 
are the eigenvalues of the local dichotomic observables, 
±1. Finally is a certain parametrization of the local 
setting di. It is introduced to facilitate further manipu- 
lations. 

On the other hand, the A^-qubit correlation functions, 
in quantum theory, have the following form: 



E{di, aat) = T • (fli ® ... ® qn), 



(2) 



where T is the correlation tensor for a quantum state, p, 
with components given by 



(3) 



with : 



Ik 



1,2,3 representing some three (local) basis 



vectors for the fcth observer. By the symbol • we repre- 
sent the scalar product in R^^ . Note that the left hand 
side of Eq. (|2Jl gives the general form of the class of cor- 
relation functions of a rotationally invariant form men- 
tioned in the introduction. If one constrains the measure- 
ment settings of each observer to just one plane (which for 
each of them can be different) the measurement direction 
vectors can be expressed by di{ai) = cosaiiji + sinaiXi, 
where sti , jji are two basis vectors of (which can be 
individually defined by each observer). In such a case. 
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the correlation function is a scalar given by: 
E{ai,...,aN) = ^ T'ii...,^ sin(ai + (ii - 1)^; 



.ijv = l,2 



X ...sin(ajv + (^AT - 1)-). 



(4) 



II. THE BELL INEQUALITY 

In order to show contrast between quantum and lo- 
cal realistic predictions for correlation function we use a 
simple geometrical fact: Two vectors h and q are equal 
if, and only if, their scalar product (h, q) is equal to the 
square of the norms of the two vectors: {h,q) = = 
Thus, {h,q} < ||(f|p implies that h ^ q. In our case 
the vectors are replaced by correlation functions, and we 
use the scalar product of a real Hilbert space of square 
integrable functions of a^, i = 1, A^, that is 

{h,q)^ / dai... daNh{ai, . . . ,aN)q{ai, . . . ,aN)- 
Jo Jo 

(5) 

Let us calculate an upper bound of the scalar product 
of Ehv and E. We get: 

iE„v,E) < / dai... / daNT\l{ai) 
Jo Jo 

X r,,...,„sin(ai + (ii - 



X sin(aAr + (lA, - 1)^) 

= ^ ] 2^ii...iiv 

zi ,. . . ,ijv — li2 

N „2 



(6) 



TT / dail{ai) sin (a; + {ii ~ 1)-). 



The sum © has 2^ elements containing all possible com- 
binations of sine and cosine functions. In Ref. 7] it was 
shown that 



{Ehv.E)<^''E, 



(7) 



If one replaces in the above inequality Eny by E one 
gets: 



{E,E) 



\E 



I dai... / daAr(Ti...i sinai... sinttAT 
Jo Jo 

• •• + 72...2 cosai... cosaAf)^. (8) 



Since J^^ da sin a = and J^^ da sin^ a — n we have 

(i?,£;) =^^Eh,...,,„=i,2?^'...„- This expression differs 
very much from the right hand side of Ineq. , and thus 
for many quantum states one can find strong violations of 
this bound. As a matter of fact these violations exceed, 
in the case of four or more qubit GHZ states, those that 
can be obtained with any standard Bell inequality (i.e., 
a two-setting per observer one), see 0- 
If 



r = 4 



-TV 



\E\ 



Er, 



> 1 



(9) 



the new Bell inequalities are violated. That is a corre- 
lation function E for which Ineq. © holds cannot have 
a local realistic model. The factor r is a good measure 
of noise robustness of the non-classical correlations in a 
given quantum state po- Take a state, given by 

p{V) = Vpo+{l-V)pno^se, 

where Pnoise — ^1 d, where in turn 1 is the unit operator, 
and d = 2^ stands for the dimension of the Hilbert space. 
Obviously, < < 1. Such a state violates the given 
Bell inequality by if and only if y > 1/r. 



III. PARTIAL SEPARABILITY 

Let us discuss how does the inequality fare in the con- 
text of discretion between partially separable states and 
states that have true N qubit entanglement. We start 
first with the case of biseparable states, for which the 
reasoning will be presented in some detail. Next we dis- 
cuss generalizations. 



A. Biseparability 



max is the maximal possible value of a correlation 



where E, 

function, that is Emax — maxo 

is important to stress that the maximalization is only 
within some fixed planes of observation defined for each 
observer by his/her local versors Xi and jji. Note that 
this inequality depends on the correlation function for 
quantum state to be analyzed. Interestingly, one can use 
also any correlation function that has the rotationally 
invariant form, as the one given in Eq. Q . This general- 
ized Bell inequality, Ineq. 0, is the basis of our further 
considerations. 



Our definition of biseparability runs as follows. As 
E{ai, . . . , un). It biseparable we shall consider all states, pj**"^, for which 
there exist a convex expansion into density matrices, pi, 
where / is some index, such that every density matrix 
Pi factorizes into at least two density matrices of two 
separate subsystems, one consisting of N — k qubits, and 
the other of K qubits, i.e. pj — pA, ® PBi (subsystems 
A J and Bj may vary, depending on /). For example 
for iV — 4 we can have the decomposition in which the 
following classes of partially factorable density matrices 
may pop up Pi^pju and pij^pu, where i, j, k,l = 1,...,4 
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enumerate the qubits. Note that density operators of the 
form pi® P2® Pz® Pa are joint members of the two classes. 

By rbisep let us denote the violation factor of Q by 
a biseparable state. Then, if the the pi...N state violates 
the inequality ^ by a factor r such that 



r > r, 



MAX 



(10) 



whe 



„MAX 



is the maximal value ofr for any biseparable 



state, then pi...N must contain the true non-classical N - 
particle correlations, i.e. cannot be biseparable. 
We shall calculate rM^i^ , but first we show that 



ri...jv = max(||£;||V4^S„,,), 
p 

that is the maximal possible violation factor of inequality 
by any state p, is equal to . 

Let us start from the case of two qubits. Our first 
task is to find the general form the two-qubit correlation 
function, for measurements in Xi — yi planes, i = 1,2. 
The correlation function is defined as 



E{ai , as) = Tr[pcr(i) (ai) cr^^) (as)]. 



(11) 



where cr^'^\ai) = coaaiax^ - 
for the operator a'^^^ ® cr'^^ 



- sinaiCTy*'. Note that 
expressed in the zi-Z2 



product basis, only the antidiagonal elements do not 
vanish. If one denotes the element of the basis by 
{|00), |01), |10), |11)}, the non-zero elements are given by 
{ad) ® a^^^kMi~k,i-i - exp(-*(-l)'(ai + {-ly-^a^)). 
Therefore, 



E{ai,a2) = 2(cos(ai +02) Repo,0;i,i 

+ cos(q!i - Q!2) Repo,i;i,o 

- sin(ai + a2) Impo,0;i,i 

- sin(ai - as) Impo,i;i,o)- 



(12) 



Using Schwarz inequality we can show that the function 
(fT^ is bounded from above by 2(|po,0;i,il + |Po,i;i,o|)- 
Please note that this value is actually achieved by E for 
"1 = (^'o,0;i,i+^'o,i;i,o)/2 anda2 = ($o,0;i,i-'^'o,i;i,o)/2, 
where $fc,i:i-fc,i-i denotes the argument of the complex 
number Pk.i-i-k.i-i- Therefore maxQ^^Q^ ^^(ai, a2) = 
2(|po,0;i,i| + |po.i;i,o|)- This result can be easily gener- 
alized for an arbitrary number of parties. For N qubits 
the general form of correlation function for measurements 
restricted to Xi — yi planes reads: 



E{ai, ...,aN) = 2 ^ (cos(ai + (-l)'"a2 + ... + (-l)'"ajv) Re po,i2,...,jjv;i,i-J2,...,i-«N 

i2,...,iN-l=0,l 

- sin (ai + (-l)''^a2 + ... + (-l)*"aAr) Im P0,^^,...,^^■l,l-^2,■■■,l-^N) (13) 



and its maximal value is given by: 
Emax = niax E{ai,...,aN) 

Qi,...,aiv 

= 2 ^ |po,fe2,...,fcjv;l,l-fc2,...,l-'c]vl-(14) 
/i:2,...,/cjv— 0,1 

Let us express the value of {E,E) = in terms 

of the parameters of the density matrix. We simply put 
the expression given in Eq. IjlSI) into the integral form of 
Since all cosine and sine functions, which appear 



in H13|) , are orthogonal to each other (with respect to the 
scalar product defined by (01 ), and all these functions 
have the same norm, equal ^{2n)^ , the final result reads 



||£;|P = 2(27r)" ^ Kfc2,...,fc«;M-^2,...,i-.«r 

fc2....,fcjV— 0.1 

(15) 

Let us calculate the maximum of ri,,,N. We get 
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ri...N - 

2(27r)^ X)fe2,...,fc„=o,l |Po,fc2,...,/cjv;i,i-fe2,...4-fcNp 

'^f.^ j.^_Q j^\po^k2,...,kN;l,l-k2,---A-kN\ 

^ 2(27r)^ Sfc2,...,fcw=o,l |Po,fc2,...,fciv;i4-fc2,...,i-fciY P ^^g^ 
~ ^^'i^J2k2,...,kN IP0,fc2,...,fcjv;l,l-fc2,...,l-fcjv=0,lP 



The inequality in is true because the module of 

any antidiagonal elements of the density matrix satisfy 
\PanUdiagonai\ < 1/2. This Can bc easily checked for pure 
states, which, of course, give the upper bound. Thus 
the relation 2\pa„udtagonai\'^ < \panUdiagonai\ uscd in the 
inequality is justified. The obtained bound is actually 
achieved by N qubit GHZ states. Simply, the correla- 
tion tensor T of a iV qubit GHZ state, ^(| + , . . . , + > 
+ | — — >), where |± > are eigenstates of <Jz, for ob- 
servation directions given for the i-th observer by Xi and 
iji, has 2^~^ nonzero components equal ±l's. Therefore 

The maximum r^f^^ can be calculated in the same 
way as above. To this end, note that the modulus of 
any antidiagonal element \p'anUdtagonai I biseparable 
density matrix is at most equal to 1/4. This again can 
be checked by considering pure biseparable states. Thus, 

now we have 4 1 Pa„tjdjagonaJp < \PanUdiagonal\- Usiug this 

fact in lfTB|) we get 



MAX _\\E'^''''''A? 1 



-— _ ^(7r/2)^. 



(17) 



Finally, the condition (|10|l one can rewritten in the form 



B. Other partial separabilities 

The presented method can be simply generalized to 
the case of fc-separable states (fc > 2). We shall use 
the following definition of A:-separable density matrices, 
^k-sep ^ There must exist for it a convex expansion into 
density matrices, p/, where / is some index, that is 

Pi.T.TV^ = J2iPiPii "^ith Pf > and J^i^Ij such that 
every density matrix pi factorizes into at least k den- 
sity matrices of k separate subsystems S^, n = 1, . . . ,k 
(such subsystems can consist out of of minimum one up 
to iV — fc + 1 qubits). That is, pi — (^n=iPs' ■ course, 
the splitting into subsystems may vary, depending on /). 
Obviously, for any antidiagonal element of fc-separable 
matrices we have \pant7lagonai \ < (1/2)^ Therefore 



„MAX _ 2- 



N 



For fully-separable states, i.e. 



(19) 



one has 



„MAX 



• fuliy-sep = i^/'^)'^ < 1- Again the inequahty 10 does 
not have to be violated to show that the state is not 
fully-separable. Please note that for all N and fc the 

^MAX 
k — sep 



threshold violation factor r^l'^^, is two times bigger than 



^MAX 
(fc— 1) — Sep* 



IV. DISCUSSION 



r(V2) 



N 



<r<-{n/2r. 



(18) 



Please note that value of upper bound is two times higher 
than the lower one. That is, the inequality (TJ, which we 
use for the construction of the condition is more 

sensitive than the ones proposed in The same sen- 
sitivity factor was obtained in however, it requires 
a specific choice of observables. Interestingly, for N = 2 
and N = 3 any violation of the inequality ((JJ) implies that 
the state has the true A''-body correlations. As a matter 
of fact, the state has true ./V-qubit correlation whenever 
r of {T)) is greater than 0.62 and 0.97 for two and three 
qubit states, respectively. That is, the inequality does 
not have to be violated. 



As we have seen the approach is pretty efficient in dis- 
covering quantum correlations of truly N particle nature. 
What is important, to get the data needed to find the 
rexp, all one needs to measure are the values of the corre- 
lation functions at two settings per each observer, exactly 
like in the case of the standard Bell inequalities. The 
optimal choice is that the two settings of each observer 
define complementary measurements. Having these it is 
a straightforward exercise to find E„iax ■ One simply uses 
the formula Q with Ti^...iM = E{xij^, ...,Xt^)exp, and 
searches for its maximal value. Please note, that exactly 
this method is implied by the derivation of the bound 
in O . Behind all this is the assumption of rotationally 
invariant form of the correlation function, which is well 
justified by the Noether Theorem, and its implications. 
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The other feature of the present approach is that it 
is very easy to formulate the general condition for an 
arbitrary number of qubits, N. Also it can be generalized 
to different type of partial separability than just the one 
involving two subsystems (bi-separability) . 

Let us finish with a purely speculative conjecture. It 
might be the case, that a conjunction of rotational invari- 
ance principle, which must hold for the form of physical 
laws, of local realism, and of some other principle, per- 
haps the one introduced by Zeilinger and Brukner 0, 



may lead to an ultimate version of Bell's theorem for 
qubit systems. One can put forward this speculation, be- 
cause more generalized inequalities involving rotational 
invariance, of the type introduced in lead to even 
more sensitive conditions for true non-classical N qubit 
correlations than which in case of some specific ex- 
amples are equivalent to non-separability criteria 1^ . 
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